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Abstract 

We study multiple scatterings of jets on constituents of quark gluon 
plasma and introduce energy-energy correlations to quantify their effects. 
The effects from a longitudinally expanding plasma on medium as well as 
high energy jets are found to be significant at both RHIC and LHC energies. 
Because jets escape from the plasma long before the completion of mixed 
phase, these effects are free from complications of final state hadronic inter- 
actions and decays. These suggest that jets can be used to probe the plasma 
that might be created in future high energy heavy ion collisions. 

PACS numbers: 25.75. +r, 12.38.Mh, 24.85. +p, 12.38.Bx 
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1 Introduction 



In high energy nuclear coUisions new states of matters, e.g. the quark gluon 
plasma predicted by numerical calculations of lattice gauge theories, might 
be produced. Because of the many-body heavy-ion system and soft multi- 
hadronic production, major efforts are to find experimental observables that 
could provide information on the states of matters produced in the collisions 
[1]. Many observables that have been intensively studied are mostly prod- 
ucts created in the plasma, like dileptons, photons, J/ip and strange particles. 
These processes can carry information about the plasma, but some of them 
are usually complicated by final state hadronic interactions and decays. 

In this paper we are interested in jets propagating through the plasma. 
It is of great interest for following reasons, (a) Jets have been extensively 
investigated both experimentally and theoretically. Because they are pro- 
duced by hard collisions of quarks and gluons, they can be and have been 
calculated with perturbative QCD to various degrees of accuracy. The agree- 
ment between theory and experiments shows that jets are one of the very 
few well understood high energy strong interacting phenomena. (b)Jets are 
produced at the very beginning of the collision process and they propagate 
in the plasma and interact strongly with the plasma constituents before they 
escape. Because of these properties they can be used to probe the plasma, 
(c) Since the final state hadronic interactions and decays take place after the 
phase transition that is expected to complete after the jets escape from the 
plasma, the jet signals are virtualy free of background of final state interac- 
tions and decays. (d)Jets are easily accessible experimentally. The typical 
jet energy is greater than 10 GeV which is much larger than the typical ener- 
gies of constituents in a thermalized plasma at the temperature that can be 
created in the future heavy ion collider, such as the RHIC and the LHC. Jets 
are therefore expected to escape from the plasma as well identified objects. 

The energy-energy correlations of hadrons produced in e^e~ collider were 
suggested by Basham, Brown, EUis and Love [2] to test QCD, and were in- 
vestigated by a number of experimental groups [3]. In hadronic and nuclear 
collisions jets are produced by hard scattering of partons. In a 2 2 process 
the two jets are produced back-to-back in their rest frame. Thus the QCD 
contribution to the energy-energy correlations is dominated by 2 — > 3 pro- 
cess. Ali, Pietarinen and Stirling [4] showed that transverse energy-energy 
correlations depend very weakly on the structure function and the normal- 
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ized correlation function is approximately proportional to coupling constant 

A jet produced at the beginning of a collision may scatter many times on 
plasma constituents before it escapes from the plasma. The average trans- 
verse momentum gained in the multiple scattering can be as high as a few 
GeV. This effect on a jet of energy of several tens GeV is of the order of 
as- In fact, as we shall see the multiple scattering has significant effects on 
energy correlations and dominates over the QCD contribution in some cases. 

The paper is organized as follows. In section 2 we start with definition of 
energy-energy correlation and derive a general formalism for evaluating the 
scattering effects on correlations. We then calculate these effects in a longitu- 
dinal hydrodynamical model that describes plasma expansion. Hadronization 
effects are discussed at the end of the section. In section 3 we present the 
numerical results and discuss their dependence on colliding energy, jet energy 
and nuclear size. In section 4 we summarize our conclusions. 

2 Jets in quark gluon plasmas 

In nuclear collisions at RHIC and LHC energies, jet productions are no longer 
rare events. High transverse energy jets stand high beyond the associated 

soft productions and can be readily identified experimentally. Low transverse 
energy jets arc, on the other hand, burried under the large amount of soft 
processes. To reduce the soft background, one studies high transverse energy 
jets. To see the effects of plasma on high Ex jets, it is essential to find 
sensitive variables. In the following, transverse energy-energy correlations 
are introduced to study such effects. As we shall see that the plasma indeed 
has significant effects on the correlations. 

2.1 Energy— energy correlations and effects of plasmas 

As aforementioned, the energy-energy correlations were introduced in [2] to 
test the QCD in e~^e~ collisions, and transverse energy-energy correlations 
were suggested in [4] to measure ag in pp collision. In nuclear collisions we 
can similarly define the transverse energy-energy correlations for high Ex 
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jets as 



^^^^^^ " § ^^^)e.ent (l) 

Here En is the transverse energy of particle i, S = Jq'^ {dTi/ d(j))d(j)^ and Acj) 
is the angle between the particle pair (i, j) on the transversal plane. The 
(• • ■)event mcans average over many events. To obtain higher statistics one 
might integrate over a range of transverse energies. In this case one has 




The numerical calculation of correlation function has been done in ref. 
[4] based on the perturbative calculations of hard 2 — > 2 and 2 — > 3 processes 
in QCD [5]. 

Now let us consider the effects of the plasma. Because two-jet pro- 
duction dominates over multiple-jet production we need only to consider 
the plasma effects on two-jet events. Let and p'^ be the four momenta 
of a jet before and after the multiple scattering off plasma constituents. 
In high energy elastic scattering, the scattering cross section depends only 
on four-momentum transfer squared, or the transverse momentum because 
— (p^ + p'fj)"^ = 4p^ sin^(r7/2) ^ p'^. The px distribution of the jet de- 
creases rapidly with increasing px and vanishes in the backward hemisphere 
{rj > 7r/2). We assume that fpiPr) is the px distribution of one jet which 
we define as the forward jet and fsiPT) is that of the other jet in backward, 
and they satisfy the following normalization condition 

/ ^iFiPT) = J ^fsiPT) = 1. (3) 

Here E — Q/2 and Q is the total jet energy. We note 

/ ^/Hpt) = I £ d'pxMpT) jf^^ dp,. (4) 

After integrating over p, one can expand E'^ — px as —p\/2E'^ -\ ). 

As we shall see later that the effects of the plasma on jets are of the order of 
{px) I E, we can neglect {px) / E"^ and higher terms. That is 

j d'^pxfpiPT) = j d^pxIsipT) = 1- (5) 
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Let da/dQo be the cross section for two-jet production in a nuclear col- 
lision, with Qq being the usual solid angle. The double energy cross section 
after the scattering on plasma constituents can be written as 



dnidn2 

which can be rewritten as 



After replacing p with pt/ sin 77 (77 is the scattering angle with respect to the 
original jet axis) in the forward hemisphere and p with px/ sin(7r — rj) in the 
backward hemisphere, we obtain 

"^'^ ^ / rf^o^ [Fivi)Fin - V,) + Fin - Vi)Fiv2)] (8) 



dflidn2 2 
where 



F{r]) = ^ Jp^dpfpiPT), for r/ < f (9) 



From fRI) we see that 



J dQFir]) = 1. (10) 



Using relation p = px/ sinr], (|^) can be written as 

E sin^ ?7 7o 

When sin 77 >> {pt)/E, we can use (|) and obtain 



^(^) = l^r-r- L PTdPTfriPT). (11) 



m = TTW^^ for < f (12) 
ZttF sm ?7 

Here (pr) is the average transverse momentum defined by 

(Pt) = J PTd^VrfM- (13) 
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The integration in can be carried out numerically [6]. To obtain a 
simple analytical result, we note that in nuclear collisions only high trans- 
verse energy jets can be well identified. Also the study of jet propagation in 
longitudinal direction is very complicated by the fact that plasma expands 
in the same direction. For these reasons we restrict ourselves to the study of 
transverse jets. Near the transverse direction the da/dQo varies slowly with 
dflo, while F[rj) varies rapidly near rj = 0. Thus the dominating contribu- 
tions to the integration in (H) come from small angular regions about the two 
detection directions Qi and As an example we consider the contribution 
from integration over the cone about Qo ~ Qi which implies rji 0, 772 ~ 4>- 
Here is the angle between Qi and ^2 on the transverse plane. Since the 
F{r]i) changes rapidly, we can expand da/dQo in the integrand as a Taylor 
series in rji for fixed 0. The dominating contribution to the integral over this 
region is then obtained by taking dQo ~ dQi, da/dQo ~ da/Qi and ri2 ~ 0. 
The integral over AQi becomes 

-^F(7r-0)/ dn,Fir^,) = -^Fin-<P), for > f (14) 

Similarly, the integration over other cones can be evaluated. Summing over 
all contributions we obtain the double energy cross section 



rf^S 1 

-F 77 - 



dnidn2 2 



da da 

+ 



dQi dQ2 



for > f (15) 



Integrating over double solid angle and keeping constant we obtain the 
energy-energy correlation 



1 2{pt) 1 



for > f , (16) 



Srf0 ttQ sin3(7r-0)' ^ 2' 

and there is no correlation for < 7r/2 in the leading order approximation. 
Here Q is the total jet energy on the transverse plane. The asymmetry about 
= 7i/2 comes from the following fact. Because the fpiPr) and fsiPr) 
decrease rapidly with increase scattering angles with respect to two initial 
back-to-back jets, the correlation is the strongest at = vr and decreases 
monotonously with decreasing and eventually vanishes at = 0. The 
discontinuity of about = 7r/2 occurs because of the approximation. To 
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avoid the discontinuity we can extrapolate the correlation to the region (j) < 
7r/2 with 

This is sensible because in < 7r/2 the correlation decreases monotonously 
to zero at = as it should be, and the correlation function as well as its 
first derivative are continuous at = 7r/2. 

2.2 Multiple scattering in expanding plasmas 

Now we turn to the calculation of (pr)- Let g{\<.T) be the normalized trans- 
verse momentum distribution of a jet after one scattering with a plasma 
constituent and P{n) be the normalized probability that a jet scattering n 
times with the constituents. The transverse momentum distribution of a jet 
after multiple scattering with plasma constituents is then 

/n n 
n d'k^i5(kTi) • • ■9{\^Tn)KVT - EkT,), (18) 
n=u j=l j=l 

with the understanding that /(pr) = (^(pt) when n = 0. In an equilibrium 
plasma the multiple scatterings are independent from one another. In this 
case one can easily show 

(pi.) = iV(4), (19) 

where (p^) is the average transverse momentum after multiple scatterings 
defined as 

(Pt) = / d'pTplfFiPT), (20) 

and {k^) is the average transverse momentum after one scattering defined as 

(4) = J d^kT^rgi^T), (21) 

and is the average number of scatterings given by = I]n=o^-P(^)- 
Noting that the (pf.) is proportional to (pr)^, and {k^) to (A;t)^, and both 
distribution should be similar in shape, we expect 

(Pt) = ^{kr), (22) 
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In hydrodynamical models can be calculated from 

(23) 

Here Ui is the particle density of the ith plasma constituent, ai is the scat- 
tering cross section of a jet with the ith constituent. The integral is carried 
out along the jet trajectory. 

We consider central nuclear collisions where the created plasma is con- 
strained in a cylinder of radius R = 1.2A^^^ (fm) where A is the nuclear 
mass number. The evolution of the plasma after its formation is described 
by hydrodynamical equations. In the Bjorken model [7] the system expands 
isentropically along the longitudinal direction and the entropy density de- 
creases according to 

sir) = s(r,)r,/r (24) 

while the temperature falls according to 

T{T)=T,[n/rf' (25) 

where is the proper time at which the plasma is formed and Tj is its 
initial temperature. When the temperature drops to Tc at time Tc a mixed 
phase emerges. In the end of the mixed phase at time Th the system is fully 
converted into hadrons. It can be shown that 

Th = TTc (26) 

where r is the ratio of degree of freedoms in the plasma phase to that in the 
hadronic phase. 

Now let us rewrite (|^) as 

N = apiJ drnir) (27) 

where api is the averaged cross section and n(r) is the total particle density 
in plasma. In the additive quark model the quark-quark crosss section is 
aqq ~ 4 mb and in perturbative QCD one has agg = (9/4)(jqg = (9/4)^(Tqg. 
Hence, we might expect api ^ {lQ + A.QNf)/{lQ + 12Nf)agg ^ 12 mb. Noting 



N 



dx 



ndx ai 
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that for an idea gas the particle density is proportional to the entropy density 
and using the expansion equation (El) we have 



(28) 



where n(rj) is the initial particle density. Tf is the time for the jets to escape 
from the plasma, if this happens before beginning of the mixed phase, i.e. 
Tf < Tc- When Tf > Tc jets will interact with the mixed phase constituents: 
quarks, gluons and hadrons. 

Because of the large degree of freedom in the plasma phase (c.f. eq. |26|) , 
the mixed phase is long and jets propagating transversal to collision direction 
are expected to escape from the plasma before the completion of mixed phase. 
Thus hadronic gas has no effects on jets. 



We note ag.^ = {2/9)app ~ 9 mb. For simplicity we assume agj, = api. 

^ . There 



In this case the number of collisions can still be expressed by (^ 
Tf is the escape time of jets which may have been scattered by mixed phase 
constituents. When a jet pair is created at a distance Ro from the collision 
axis, the average number of collisions each jet experienced is 



N = -(ypinn^Ti) In 



it — Jig 



(29) 



and iV = if 



Rl 



< Tf. Since Rq is unmeasurable, an average over 



geometry must be made. For a sharp nuclear density profile we obtain 



Noting that n{Ti 



N = ^apiTin{Ti) 



s{Ti)/c with c 

,2 



R' I t: 
In ^ - - + 



T 



(30) 



3.6, we obtain 



4:5c 



In- 



R^ 



1 tI_ 

2 ^ 2i?4 



(31) 



where qq = IQ + 21Nf/2 is the degree of freedom in plasma and Nf is the 
number of active quark flavor. 

In the above calculations only longitudinal expansion is taken into ac- 
count. When the system expands transversally as well, it will take longer for 
jets to escape from the plasma. This may increase the number of collision and 
enhance plasma effects. On the other hand, transverse expansion accelerates 
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the cooling process and reduces plasma density more rapidly. Hence, the re- 
sulting effects is expected to be small. Furthermore, we note that transverse 
expansion is negligible in the early stage of the expansion. As a result, the 
transverse expansion has little effects on jets. 



2.3 Hadronization 



The results obtained so far are for quarks and gluons. To compare with exper- 
iments, the effects of hadronization should be taken into account. Basham et 
al. [2] have calculated that the hadronization correction to the energy pattern 
and correlations. Similarly one can show that the hadronization correction 
to the transverse energy-energy correlations is given by 



Co(M 
ttQ sin'^ ( 



for 6 < TT. 



(32) 



Here cq = 2.5 is the mutiplicity density at rapidity y = and (/it) is the 
average transverse momentum of hadrons from jet fragmentation [2]. Because 
two hadrons with angle (p may come from either the same jet or from the 
two separate back-to-back jets, the hadronization correction to correlation is 
symmetric about = 7r/2. 

Summing over the effects of QCD corrections, multiple scattering in 
plasma and hadronization, we finally obtain the energy-energy correlation 



1 rfE(«^^) 



f 2^{kT) + Co(/ir) 
TxQ sin^(7r — 0) 



ttQ 



+ 



nQ sin^ < 



for (f) > 7r/2 
for < 7r/2. 



(33) 



3 Results and discussion 

From (^) and (|3TD we see that the results depend on three parameters, the 
initial temperature Tj, the jet energy Q and the geometrical size R, that can 
be changed in experiments. 

We first look at central Au-Au collisions at RHIC energy (200 GeV per 
nucleon) where the initial temperature is expected to be about 250 MeV. The 
results for correlations are shown in Figs. 1 at jet energies Q = 10, and 20 
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GeV. In the calculation we used Nf = 3, api = 12 mb, A = 0.2 GeV, = 1 
fm, {kr) = {hr) = 0.36 GeV and the results of refs. [4, 5] for dE^^^^^dcf). 
As we can see from the figure that the hadronization effects are unimportant 
at Q = 20 GeV, while the plasma effects are distinctly beyond the QCD 
effects. To minimized the hadronization effects one studies the asymmetric 
correlations defined as [2] 

for < 7r/2 (34) 

In Fig. 2 we show the corresponding A{(p) distributions. Prom the figure we 
see that the plasma effects are significant and should be easily measured ex- 
perimentally if the plasma is created in the collision. Here the hadronization 
effects cancel because of the symmetry about = 7r/2. 

In Figs. 3 and 4 we show the {l/T,){dT,/d4>) and A{(j)), respectively, for 
Q — 40, and 100 GeV at LHC energy (6300 GeV per nucleon) where we esti- 
mate Tj ~ 500 MeV. Compare Fig. 1 with 3 and 2 with 4 we see that plasma 
effects increase with increasing initial temperature. At LHC the plasma ef- 
fects stand well beyond the effects of QCD corrections and hadronizations 
for jet energy as high as 100 GeV. 

To see the plasma effects more directly, one might compare the integra- 
tions that are defined as 

1 dT, 

'^-L '''' 

for energy-energy correlations and, 

a A = 1'^' A{4>)d(t> (36) 

for asymmetric energy-energy correlations. The results for various Tj and Q 
are listed in Tables 1 and 2 respectively. Here we integrated over (f) from 7r/6 
to 57r/6 for cxc, and from 7r/6 to 7r/2 for cta. From the tables we see that the 
integrations are enhanced by a factor of 1.2 to 1.6 due to the presence of the 
plasma. In e^e" experiments the integration can be very well determined [8]. 
Similarly, we expect that such enhancement can be unambiguously measured 
in high Q jet events in heavy ion collisions. 
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Table 1: Integrations of energy-energy correlation function for various Tj and 

Q 



T\ (MeV) 


Q (GeV) 


QCD 


QCD+Hadr. 


QCD+Hadr. +Plasiiia 


250 


10 


1.29 


1.54 


2.27 


250 


20 


1.32 


1.44 


1.80 


500 


40 


0.93 


0.99 


1.51 


500 


100 


0.95 


0.98 


1.19 



Table 2: Integrations of asymmetric energy-energy correlation function for 
various Tj and Q 



Ti (MeV) 


Q (GcV) 


QCD 


QCD+Hadr. 


QCD+Hadr. +Plasma 


250 


10 


0.76 


0.76 


1.12 


250 


20 


0.33 


0.33 


0.52 


500 


40 


0.54 


0.54 


0.80 


500 


100 


0.24 


0.24 


0.35 



In central collisions the initial particle density increases with A^^^. When 
R » Ti the number of scatterings increases Wiili A^/^\n{R^/Tf). In Figs. 
5 and 6 the (l/S)((iS/d0) and distributions for Q = 40 GeV for S-S 
collision are compared with those for Au-Au collision at LHC, respectively. 
We see that the results depend rather weekly on the nuclear geometry. 

We note that in the kinematic region that we are interested in, domi- 
nating effects arc from scatterings of two gluon jets, and our results are not 
applicable to scattering angle sin?7<(pr)/Q which corresponds to the region 
of and (f) ^ n. 

In small scattering angle the multiple soft gluon bremsstrahlung plays an 
important role [9,10]. In the regions of ~ and (p ^ n one might study 
the momentum imbalance or acoplanarity as suggested by Appel [9] , Blaizot 
and McLerran [10]. 
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4 Conclusions 



We used (asymmetric) transverse energy-energy correlations to study the 
effects of multiple scatterings of jets on constituents of quark gluon plasma. 
In the calculations we considered the expansion of plasma, QCD corrections 
and nonperturbative effects. Our results show that the scattering effects 
are remarkably strong in (asymmetric) correlations. At LHC they stand 
well above the background from QCD corrections and hadronizations for jet 
energy as high as 100 GeV. If the plasma is formed in nuclear collisions and 
go through a phase transition, the jets are expected to escape before the 
completion of mixed phase due to the large degree of freedom in plasma. 
In this case the observed effects are free of final state hadronic interactions 
and decays. These imply that jets could be used to probe the quark gluon 
plasma. 

The plasma effects increases with decreasing jet energy. To obtain maxi- 
mal effects, jet energy should be kept as low as possible. We note, however, 
that low energy jets, or minijets, can be created or quenched due to mul- 
tiple interactions with plasma constituents. Because of fluctuation in soft 
processes [11], it is also elusive to identify minijets in experiments. To avoid 
such backgrounds and experimental ambiguities, jet energy should be kept 
high enough. Total jet energies of Q = 10 and 40 GeV might be s ensible 
thresholds at RHIC and LHC energies, respectively. 

The main background of jet signature is probably the hadronic gas. When 
jets propagate in hadronic gas, the multiple scattering between jets and 
hadrons results in similar effects as plasmas do. In general one expects that 
hadronic gas exists at lower temperature. In very high energy heavy ion 
collisions, such as RHIC and LHC, the effects from hadronic gas should be 
smaller than those from plasmas. To make quantitative comparison, realistic 
models and detailed calculations are necessary [6]. 
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Figure captions 

Fig. 1 The (l/S)(c^E/c^(/)) distribution for Tj = 250 MeV in Au-Au colli- 
sions, (a) for total transverse jet energy Q — 10 GeV, and (b) for 
g = 20 GeV. 

Fig. 2 The A{(f)) distribution for the same parameters as those in Fig. 1. 
(a) for Q = 10 GeV, and (b) for Q = 20 GeV. 

Fig. 3 The (l/E)(dE/d0) distribution for Ti = 500 MeV in Au-Au colli- 
sions, (a) for Q = 40 GeV, and (b) for Q = 100 GeV. 

Fig. 4 The A{(f>) distribution for the same parameters as those in in Fig. 3. 
(a) for g = 20 GeV, and (b) for Q = 100 GeV. 

Fig. 5 The comparison of (l/S)((iS/(i0) distribution for S-S and Au-Au 
collisions for Q = 40 GeV at Ti = 500 MeV in Au-Au coUisions. 

Fig. 6 The comparison of A{(j)) distribution for S-S and Au-Au collisions 
for the same parameters as those in Fig. 5. 
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